On the nature of highly vibrationally excited states of Thiophosgene 



Srihari Keshavamurthy 
Department of Chemistry, Indian Institute of Technology, Kanpur, Uttar Pradesh 208016, India 

In this work an analysis of the highly vibrationally excited states of thiophosgene (SCCI2) is made 
in order to gain insights into some of the experimental observations and spectral features. The states 
analyzed herein lie in a spectrally complex region where strong mode mixings are expected due to 
the overlap of several strong anharmonic Fermi resonances. Two recent techniques, a semiclassical 
angle space representation of the eigenstates and the parametric variation of the eigenvalues (level- 
velocities) are used to identify eigenstate sequences exhibiting common localization characteristics. 
Preliminary results on the influence of highly excited out-of-plane bending modes on the nature of 
the eigenstates suggest a possible bifurcation in the system. 



I. INTRODUCION 

Understanding the nature of the highly excited molec- 
ular eigenstates is equivalent to deciphering the mecha- 
nism of intramolecular vibrational energy redistribution 
(IVR) occuring in the molecule [T]. However, the assign- 
ment of eigenstates is far from simple. The existence 
of and interplay of several strong anharmonic resonances 
result in complicated spectral patterns and highly mixed 
states. Nevertheless, eigenstates of molecular systems 
rarely exhibit the extreme scenario of being completely 
ergodic or completely regular [5]. A generic situation, 
even at fairly high energies, is the coexistence of sev- 
eral classes of eigenstates with differing degree of mixing 
being interspersed among each other |3ii4j. More interest- 
ingly, the interspersed states can show several sequences 
associated with certain identifiable common localization 
characteristics. In other words, although a full set of 
quantum numbers do not exist for labeling such states, 
there are reasons to believe that a sufficient number of 
approximate quantum numbers do exist to understand, 
and hence assign, the eigenstates. Since the approximate 
quantum numbers arise out of the local dynamics due 
to specific resonances, relevant at the energies of inter- 
est, the assignment is inherently dynamical in nature. 
In a nutshell, several decades of work have shown that 
dynamical assignments can be done reliably only if the 
structure of the corresponding classical phase space is 
well understood 0. 

This work is concerned with the dynamical assignment 
of the eigenstates of thiophosgene SCCI2. The Hamilto- 
nian of interest is a highly accurate spectroscopic Hamil- 
tonian obtained[6] by Gruebele and Sibert via a canoni- 
cal Van Vleck perturbation analysis of the experimentally 
derived normal mode potential surface[7]. The Hamilto- 
nian can be expressed as H = Hq + T^es with 




Vrcs = K<i26a2a5al + Ki<ieaialal 



+ K23ia\a2al + K26ia{a2al + c.c. (1) 

The zeroth-order anharmonic Hamiltonian is Hq, 
whose eigenstates are the feature states and K-os accounts 
for the important anharmonic resonances responsible for 
the IVR dynamics. In the rest of the paper the vari- 
ous resonances will be refered to by simply indicating 
the modes. For instance, the first resonance term in V^os 
above will be refered to as the 526-resonance. The hermi- 
tian adjoints of the operators are denoted by c.c, and the 
operators ak, a\, and Vk = ajjflfc are the lowering, raising, 
and number operators for the fc**^ mode respectively with 
[afc,a^] = 1, in analogy to the usual harmonic oscillator 
operators. Based on the possible resonances suggested 
by the normal mode frequencies and the observed spec- 
tral "clumps" in the dispersed fiuorescence spectra, three 
approximately conserved quantities (known as polyads) 
arise which can be expressed as 

K ^ Vi + V2 + V3, 

L ^2vi+V3 + V5 + vq, 
M = W4. (2) 

Thus, the i>4 mode is not involved in any resonant inter- 
actions with the rest of the modes. Due to the polyad 
constraints one has effectively a three degree of freedom 
system and three of the quantum numbers out of a total 
six needed to assign a quantum state are already known. 
Nevertheless, as will be seen below, even this reduced 
system presents a stiff challenge with regards to the dy- 
namical assignments of the highly excited eigenstates. 

The IVR dynamics in SCCI2 has been the subject of 
several recent experiment al[7HS] and theoretical (TDHH] 
studies. In an early detailed study, Strickler and Grue- 
bele estabhshed[7] that the IVR state space (zeroth-order 
quantum number space) has an effective dimensional- 
ity of three as opposed to the theoretically possible six. 
Combined with the observation of regular progressions 
persisting upto the first dissociation energy of about 600 
THz, one expects restricted IVR in the system. This 
in itself is surprising since the effective Hamiltonian fit, 
as in eq. [l] to the stimulated emission pumping (SEP) 
spectra contains strong multimode Fermi resonances. In 
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addition, the SEP spectra are quite complex and con- 
gested in the 200—300 THz region. However, it was noted 
that the regular progressions mainly involved the Franck- 
Condon active i^i (C-S stretch) and the (out-of-plane 
bend) modes, and thus avoided the strongest resonances. 
Furthermore, Chowdary and Gruebele recently[8] showed 
the existence of sharp SEP spectral features in the en- 
ergy range of 550 — 700 THz which includes both the 
molecular SCCI2 ^- CS -f CI2 and the radical SCCI2 -> 
SCCl -I- CI dissociation channels. Although, one could 
rationalize the observation as due to the avoidance of 
the strong 526-resonance involving modes 1^2 (ai symme- 
try C-Cl stretch), (62 symmetry C-Cl stretch), and 
(^2 symmetry Cl-C-Cl bend), a similar argument for 
the equally strong 156-resonance involving modes vi (ai 
symmetry C=S stretch), 1^5, and vq is not possible. Sib- 
ert and Gruebele performed :6| a detailed study of the 
IVR dynamics using eq. [T] and concluded that around 
240 THz (- 8000 cm^^), close to the threshold for facile 
IVR, zeroth-order states coexist exhibiting a broad range 
of dynamical behvior ranging from fairly restricted to 
highly facile IVR. A rather extensive study[Tn| of nearly 
10® zeroth-order bright ("feature") states by Chowdary 
and Gruebele concluded that about one in thousand fea- 
ture states are localized even at the highest energy - not 
enough to invalidate statistical theories but certainly in- 
dicating that even near the dissociation limit SCCI2 is 
not just a bag-of-atoms. 

The above observations hint at the existence of dif- 
ferent classes of eigenstates in a specific energy range. 
For instance, even around an energy of about 590 THz 
one observes P] states whose IVR line widths vary non- 
monotonically with energy. Thus, in order to gain fur- 
ther insights into the IVR dynamics a closer look at the 
eigenstates of eq. [T]is necessary. An earlier work by Jung, 
Taylor, and Sibert using a second order version of eq. [l] 
established JJj that several distinct eigenstate classes ex- 
ist in the spectrally complicated region 7000 — 9000 cm~^ 
(- 210 - 270 THz). More recently [1 2 , Jung and Tay- 
lor dynamically assigned the states and came up with a 
detailed semiclassical classification scheme for the eigen- 
states. 



A. Classical limit Hamiltonian 

In order to dynamically assign the eigenstates it is im- 
portant to analyze the classical limit Hamiltonian corre- 
sponding to eq. [l]and the procedure is well established in 
literature. The full classical Hamiltonian, in terms of the 
zeroth-order action-angle variables (1,0), can be found 
in an earlier workfTT]. Here, the reduced three degree of 
freedom classical Hamiltonian, taking into cognizance the 
polyads {K, L, M), is of interest and can be obtained via 
a suitable (1,9) — ^ {J,ip) canonical transformation[TT] . 
The resulting reduced Hamiltonian can be expressed as 
H{J,iP) - i7o(J) + Kcs(J,i/') where 

3 3 

dij Ji Jj 

i—1 — 1 

3 

+ ^ PijkJiJjJk, (3) 
i^j.k—1 

is the zeroth-order Hamiltonian and the resonant pertur- 
bations (nonlinear resonances) are given by 

Vros(J,'0) = vizq{3;Kc,Lc)cos.iI:i 

+ U526 (J; i^c,ic) cos ?/'2 

+ fi25(J;i^c)cos(V'i + ^2) 
+ V3e{J;Kc,Lc)cos2tp3 

+ W23l(J) COs('i/'l - -02 - 27/13) 

+ v26i{3;Kc,Lc)cos{^pi - 'ip2)- (4) 

In the above equations the actions conjugate to the an- 
gles i^Pi, V-s, V-e) are {K^ = K + 3/2, = L + 5/2, M, = 
M + 1/2) i.e., the classical analogs of the quantum 
polyads. The transformed actions {Jk} are related to the 
zeroth-order quantum numbers as Jk — (w/c + l/2)/i. The 
various parameters {C,zu,a, (3) of the reduced Hamil- 
tonian are determined in terms of the original param- 
eters (see supplementary information in reference for 
the values). The relations are not shown here but note 
that at this order the parameters (0,^7,0.) depend on 
the polyads K, L, and M. 



In this work an attempt is made to answer the fol- 
lowing questions. Firstly, can one identify the various 
eigenstate classes in a manner that is not dependent on 
visualizing the high dimensional wavefunctions quantum 
mechnically or semiclassically? More importantly, such 
a method should be able to make contact with the rele- 
vant structures like periodic orbits and their bifurcations 
in the high dimensional classical phase space. Secondly, 
do the eigenstate classes and sequences identified previ- 
ously persist when using a higher order and more accu- 
rate Hamiltonian? Finally, is there a possibility of new 
modes being born due to specific bifurcations in the sys- 
tem? 



II. DYNAMICAL ASSIGNMENT 

In this section two recent techniques O E], used to 
analyze the highly excited eigenstates, are briefly men- 
tioned. The reader is referred to an earlier publication^l^ 
for the theoretical basis and validity of the methods along 
with a detailed comparison. 

In the first approach [TTHI3] the quantum states \a) for 
a given {K, L, M) are expressed in the semiclassical angle 
space ip ^ {ilji,ilj2,^3) as 

{xP\a)= ^-e'-"^, (5) 

ve{K,L,M) 
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with V = (ui, V2, V3) and ignoring an overall phase factor. 
Due to the 27r-periodicity of the angles, the semiclassical 
wavefunctions a{i()) now "live" on the classical reduced 
configuration space which is a three dimensional torus 
T^. The nature of the states is determined by plotting 
the density Ki/jja)!^ and the phase of eq. [sjon select two- 
dimensional sections of T^. 

The second approach |141 involves computing the 
eigenstate expectation values of the various anharmonic 
resonances. The expectation values are naturally related 
to the parametric variation of the corresponding eigen- 
values, also called as the level-velocities. For example, 

BE 

{a\V^M = — ^. (6) 

Previous works [I4l - [T7] have shown that the level- 
velocities, semiclassically, are extremely sensitive to the 
underlying structures in the classical phase space. There 
is also an interesting connection (not elaborated here) 
between the parametric variation of eigenvalues to the 
matrix fluctuation-dissipation theorem [15] of Gruebele. 
In what follows the terms expectation values and level- 
velocities will be used interchangeably. 

Recently [16] it was argued that a combination of the 
above approaches can be very powerful towards gaining 
key insights into the nature of the highly excited eigen- 
states. The current work provides further support for 
such an argument. To begin with, a single resonance 
case is analyzed in order to set the stage for looking at 
select eigenstates of the full system. Although the single 
resonance case is classically integrable, the correct as- 
signment of the states needs to be done with care due to 
the multimode nature of the resonance. 



A. Integrable multimode resonance 

In the full quantum Hamiltonian, all the anharmonic 
resonances as shown in in eq. [ij are turned off except 
for the multimode 156-resonance. Since, the quantum 
and classical Hamiltonians couple only the i^i , 1^5 , and 
vq modes, the quantum numbers {v2,V3) are good quan- 
tum numbers. Consequently, one already has five good 
quantum numbers {K, L, M,V2,V3) to label those quan- 
tum states that are under the influence of the said reso- 
nance. States that are not affected by the resonance are 
trivially assigned by the full set of zeroth-order quantum 
numbers and hence will not be discussed further. In the 
present case it is desirable to organize the states using 
Ps = {L — V3) since this represents the plane on which 
the states live in the original six-dimensional state space. 

The geometry of the (fi, t^s, vq) space is shown in Fig.[l] 
for the specific case of V3 = 0. For given {K,L,AI), the 
space is filled with Pg-planes with V3 being fixed on each 
plane. On a given Ps-plane, fixed values of V2 correspond 
to lines determined by points (wi, K2 — vi,Ps — if 2 — vi) 
with K2 = K — V2- An eigenstate will be delocalized only 
along one of the lines on the Pg-plane. Two points are to 



be noted at this stage. Firstly, the area of the P^-planes 
decrease with increasing V3 and the length of the lines on 
a given plane decreases with increasing V2- Consequently, 
unperturbed normal modes should be expected for high 
excitations in the (w2, U3) modes. Secondly, the number of 
state space points (ui, K2 — vi, Ps ~ K2 — wi) on a given 
line determines the quantum number /i, an excitation 
index, related to the phase space area enclosed by the 
156-resonance island in the classical phase space. 




1 2 3 4 5 6 7 
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FIG. 1. Projection of the Ps = 14 plane onto the {vz,vq) 
space. Mixing of the zeroth-order states occurs only along 
lines corresponding to constant values of V2 (indicated) lead- 
ing to the sequences observed in Fig.|2|A) and (B). Classical 
prediction (cf. Eq. [s] grey squares) of the maximal level- 
velocities corresponding to large quantum eigenstate expec- 
tation value (Vise) of the 156-resonance operator. 

Therefore, a knowledge of the quantum number pro- 
vides a complete assignment \K, L, M, Ps]V2, fJ.) of the 
single resonance eigenstates. This information can be ob- 
tained both from the semiclassical wavefunction and the 
level- velocity approaches mentioned above. For instance, 
looking at the phase of the semiclassical state in the V'l = 
plane will immediately reveal the (u2,W3) quantum 
numbers. Similarly, keeping track of the phase advance 
along the -01 direction in the ip^ — ot 'ip2 = plane will 
yield fi as the excitation quantum number |121 113] . 

Equivalently, in Fig. [2](A) the quantum eigenstate ex- 
pectation values of the 156-resonance operator are shown. 
The apparent congested nature of the figure is a direct 
result of the interleaving of the various state sequences 
with different W3. The sequences, however, can be easily 
identified by analyzing the corresponding classical phase 
space function Ui56(J; K^., L^) cos "01 with 

Vi56^'2 [Jl {K2c - Jl ) [Psc -K2c-Jl]^'^. (7) 

In the above expression the quantities K2C and Psc corre- 
spond to the classical analogs of the quantum K2 and Ps 
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FIG. 2. Single 156-resonance Hamiltonian. (A) Quantum expectation values (Vise) as a function of the eigenenergies. The 
apparent complexity is due to the interspersing of several easily identifiable regular sequence of states. (B) (Vise) for states 
with 113 = 0, i.e., Pise = 2vi +1)5 +W6 = 14 (also shown in (A) as grey circles). Various sequences with differing V2 are indicated 
by dashed lines. A sequence corresponding to say V2 = Q (filled circles) is assigned as |7, 14, 0, 14; 0, n) with /i = 0, 1, . . . , 7 being 
an excitation index. See text for explanations. 



respectively. It can be shown that classical expectation 
value is large and has equal magnitudes at the resonant 
angles "01 = 0, tt and resonant action 



Ji 



(8) 



The expectation values are positive for t/'i = and nega- 
tive for = TT. From the classical phase space perspec- 
tive the opposite signs translate to states being localized 
about the stable or unstable fixed points of the resonance 
island. Using the classical estimate it is possible to iden- 
tify the highly localized states which further spawn the 
sequences with different ^ and in Fig. [2]jB) some of the 
?;2-sequences for fixed = Q i.e., Pg = 14 are shown. 
For instance, the sequence shown in Fig. [2]^B) (circles) 
is assigned as |7, 14, 0, 14; 0, /i) with ji — (maximum 
positive expectation value) to fi — 7 (maximum nega- 
tive expectation value) . All such sequences can be easily 
identified and hence a complete assignment of the inte- 
grable, single resonance cases can be provided in terms 
of the appropriate quantum numbers. 

Including other independent resonances in eq. [T] re- 
sults in a non-integrable system with highly mixed states 
which are presumably interspersed with some of the reg- 
ular states. For example, including the 526-resonance 
(and/or the 125-resonance) in addition to the 156- 
resonance destroys the V2 quantum number. In terms 
of the geometry described above, since V3 is still con- 



served, the Ps-planes are preserved but the states now 
spread over the entire plane. Further inclusion of the 
231-resonance and 36-resonance destroys the P^-planes 
as well and the states now spread over the entire state 
space. In the next two sections the goodness of V2 and 
W3 are broken sequentially in order to check if any new 
state sequences appear despite the nonintegrability of the 
system. 



B. Breaking the goodness of V2 

In this section only the resonances involving modes 
1^1,1^2,1^5, and are considered. Thus, V2 is no longer a 
good quantum number whereas W3 is still a good quantum 
number. Note that the set of resonances involved here are 
the strongest ones in thiophosgene and therefore strong 
state mixings are expected. Are there any patterns in 
this case? In analogy to Fig. [2](B), the W3 = case of 
the nonintegrable system is shown in Fig. [3] using the 
level velocity approach. Since several level velocities are 
compared in the figure, all of them are scaled to zero 
mean and unit variance. 



From Fig.[3]it is clear that eigenstate sequences do exist 
in this non-integrable case and a few of them are high- 
lighted by dashed lines. It is important to note that the 
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FIG. 3. Scaled level velocities for the noninetgrable system 
with 156 (thick black lines), 526 (thin lines with open circles), 
125 (black lines in top inset), and the 261 (grey lines in top in- 
set) anharmonic resonances present. Eigenstates belonging to 
the good quantum number — are shown. Two eigenstate 
sequences (dashed lines) starting at the high energy end are 
easily identified. Similar sequences at the low energy end are 
also shown. The sequences merge towards the middle result- 
ing in mixed states. (Bottom inset) Participation ratios Np 
(cf. eq. [9| in the 156 (grey circles), 526 (black open circles), 
and 125 (black line) single resonance basis. 



observed sequences are different from the D2-sequences 
seen in Fig. [2]jB) since V2 is no longer a good quantum 
number. Three key observations are worth making at 
this stage. Firstly, the initial sequences starting at high 
energy can be identified using (Vise) or equally well us- 
ing (V526), two of the strongest resonances in the system. 
Interestingly, as evident from Fig. pi top inset), these se- 
quences can also be identified from (V125) and (V261) val- 
ues. The large values of (^261) despite the weakness of the 
261-resonance is an example of the nontrivial nature of 
such multiresonant systems. The two major resonances 
induce a relatively strong 261-resonance in certain re- 
gions of the Ps plane, resulting in significant perturba- 
tion of the state sequences. Thus, at high energies the 
states are under the influence of several resonances. Sec- 
ondly, at low energies the sequences are less robust (note 
the curvature of the dashed lines in the figure) and fewer 
in number, but seem to be organized more by the 526- 
resonance. Thirdly, the various sequences "collide" near 
E ~ 7800 cm~^ leading to significant disruption of the 
regularity of the sequences. Small level velocities and 
energy spacings in this region imply several multistate 
avoided crossings and suggest that the state mixing could 



be due to dynamical tunneling [T5]. Dynamical assign- 
ments of the states in this complicated transition region 
are difficult if not impossible. 

Support for the level velocity predictions comes from 
computing the participation ratio Np in the various sin- 
gle resonance basis |r) (eigenstates of single resonance 
Hamiltonians), shown in Fig. [sj^bottom inset). The num- 
ber of basis states participating in a given eigenstate la) 
is computed as 



(9) 



The participation ratios indicate that, except for a few 
states at the high and low energy ends, most states are 
moderately mixed. Moreover, large values of Np at 
E ^ 7800 cm""'^ in every basis agrees well with the level 
velocity data. Note that hints for the existence of eigen- 
state sequences are present in the Np data, but nowhere 
as emphatically as refiected in the level velocity plots. 

The various sequences observed here are consistent [H] 
with the previous results of Jung and Taylor. This 
can be established by looking at appropriate slices of 
the semiclassical three-dimensional angle space represen- 
tations of the eiegnstates. In Fig. |4] the semiclassical 
('0ij'02;V'3 = 0) angle space representations are shown 
for select states in the V3 = manifold. A striking cor- 
respondence can be seen between the sequences revealed 
by the level velocities shown in Fig. [3] and the eigenstate 
density patterns as observed in the semiclassical angle 
space. The first column of states correspond to the se- 
quence, indicated in Fig. |3j starting with the highest en- 
ergy state and have clear nodal structure - no nodes in 
the diagonal direction and increasing number of nodes in 
the antidiagonal direction. Similarly, states in the sec- 
ond column belong to the second high energy sequence 
indicated in Fig. |3] and exhibit an additional node in the 
diagonal direction. Both these sequences correspond to 
the class A classification given earlier 12] and the nodes 
along the diagonal and antidiagonal directions yield two 
additional quantum numbers td and ta respectively. For 
instance, the first state in column one of Fig.|4]is assigned 
as \K, L, M;v3 = 0,td ~ 0,ta = 0) and the assignment 
\K, L, M; 0, 1, 2) is appropriate for last state in the second 
column. 

On the other hand, the third column in Fig. [4] corre- 
sponds well to the sequence indicated in Fig. |3] starting 
with the lowest energy state. Although the angle space 
densities show localization around ipi = 0, 27r, suggest- 
ing 156-resonant states, the level velocity results show 
that the infiuence of the 526-resonance is far greater in 
comparison. State space plots and the participation ra- 
tio data do establish the first state of the sequence to 
be closer to being a 526-resonant state. However, the 
later states in the sequence are perturbed strongly by the 
induced 261 resonance (see the top inset of Fig. 3 ) and 
hence nontrivial to assign. The last column of Fig. 4 show 
three states that are in the complicated energy region 
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FIG. 4. Semiclassical (i/;!, ■02) angle space representations of select states corresponding to the case shown in Fig.jsjwith V3 = 
being a good quantum number. The "03 = slice of the full three dimensional space are shown. In every plot the angles range 
from to 27r. (First column) States belonging to the first sequence in Fig. [s] starting at the highest energy. (Second column) 
States forming the second sequence in Fig. |3] at high energy. (Third column) States belonging to the first sequence in Fig. [s] 
starting at the lowest energy. (Last column) States that do not form any clear sequence. See text for discussions. 



around 7800 cm^^. The mixed nature of these states is 
evident both from the angle space densities as well as the 
level velocity and participation ratio data. Presumably 
the induced 261-resonance is leading to the highly mixed 
states. In particular, the middle state of the last column 
seems to be a state that is delocalized over an extended 
region of the phase space and such states have been ob- 
served in a different context in an earlier studvfTB]. 



C. Breaking the goodness of V3 - the full system 

The previous section identified sequences of eigenstates 
exhibiting similar localization properties of a subsystem 
of the full Hamiltonian eq.[l] wherein was a good quan- 
tum number. The goodness of ^3 is strictly broken in the 
full system, obtained from the subsystem with the inclu- 
sion of the remaining 36 and the 231 resonances. A key 
question is whether the sequences in Fig. [3] would persist 
in the full system. The answer to this question can be ob- 
tained by analyzing Fig. [s] which shows scaled (Vise) and 
(V526) for all the 288 eigenstates. The expectation values 
(V261) and (V125) are also substantial, as in Fig. [s] but 
not shown for the purpose of clarity. It is immediately 
clear from Fig. [5] that the high end of the polyad does ex- 
hibit several sequences highlighted by solid and dashed 
lines connecting the states. Specifically, the solid line 
sequences correspond to increasing W3 with fixed {td,ta) 



quantum numbers and the dashed lines correspond to in- 
creasing tj^ with ta — and fixed W3 values. In particular, 
these high energy states correspond precisely to the class 
A states discussed in the previous section with similar 
semiclassical (V'i,'02) angle space representations. Note 
that in Fig. [3] only the U3 = states are shown and in- 
cluding states for different values of U3 would show high 
energy patterns similar to those seen in Fig. [5] 

The robustness of the class A sequences is also con- 
firmed in Fig. [sf^I) wherein the participation ratio of the 
high energy eigenstates is seen to be nearly unity in the 
W3-good basis i.e., in the eigenbasis of the subsytem dis- 
cussed in the previous section. For these set of states V3 
is approximately conserved and can be used to label the 
eigenstates. Consequently, it is possible to infer that the 
lowest frequency mode 1^3 (oi, Cl-C-Cl bend) is decou- 
pled from the rest of the modes in the [7800, 8100] cm^^ 
region. The IVR dynamics ensuing from a feature (an- 
harmonic zeroth order) state comprising mainly of the 
class A states is therefore expected to be nonstatistical. 
A crucial observation in this context is summarized in 
Fig. [sjjil) - participation ratios in the zeroth-order basis 
are unable to identify the various sequences as refiected 
in the expectation value patterns. Stated differently, de- 
localization in state space does not immediately imply 
delocalization in the semiclassical angle space and hence 
in the full phase space. 

Furthermore, two important observations can be made 
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FIG. 5. Eigenstate sequences as revealed by (Vise) (in grey) and (V526) (circles) forthe full Hamiltonian. The various classes 
of eigenstates are denoted by A,B,C,D, and E in accordance with an earlier study[T3. Sequences with increasing vs with fixed 
(td,ta) (black lines). Sequences with increasing td with ta = Q and fixed us (dashed lines). States with varying «3,td and fixed 
ta = 1 (squares). (I) Participation ratio of the eigenstates in the Vs-good basis. Note the significant breakdown of the goodness 
of in the middle of the polyad. (II) Participation ratio in the zeroth order basis of the {td,ta) ~ (0,0) sequence (solid line) 
and V3 = sequence (dashed lines). 



in the mid to low energy regions of the polyad. Firstly, for 
E < 7800 cm~^ it is clear from Fig.jsjl) that the goodness 
of V3 is broken substantially. Indeed, the regular lattice 
formed by the various sequences at the high energy end 
of the polyad begins to "dissolve" around E ^ 7800 cm^^ 
due to the influence of several resonances. Secondly, at 
the lowest energy end, the states are mainly influenced 
by the 526-resonance and hence, as noted earlier, typical 
class C states. The middle region of the polyad span- 
ning [7350,7750] cm~^ is particularly complicated with 
several highly mixed states and states exhibiting transi- 
tional nature between the two extreme classes, A and C, 
of localized states. 



III. HIGHLY EXCITED ^4 STATES 

Up until now the focus has been exclusively on the 
= case i.e., no excitations in the out-of-plane 
bending mode. Experimentally, however, the vi and 
Pi modes have dominant Franck-Condon activity[7] and 
sharp assignable features in the spectra are seen[5] even 
around and above the dissociation energy. Most of the 
sharp features have fairly large values of V4 in their as- 
signments. Thus, it is interesting to ask if the nature of 
the eigenstates, in terms of localization as well as the dif- 
ferent classes, undergo any significant change upon excit- 



ing the mode. Earlier works have hinted at a smooth 
morphing of the reduced dimensional wavefunctions0 
with increasing 1/4 excitation and the existence of IVR- 
protected zeroth-order states for large number of quanta 
in the mode [6 . Can the level- velocity technique used 
in this work provide any insights? 

As a preliminary effort. Fig. |6] summarizes the result of 
increasing 1/4 excitation on the highest energy eigenstate 
in the polyads [K — 7,L — 14, V4]. The total number of 
eigenstates remains the same (288) due to the decoupling 
of the i>4 mode from the other modes, but the eigenstates 
span higher energy ranges. The main part of Fig.[6]shows 
the participation ratios of the eigenstate in various basis. 
Clearly, the highest energy state changes it nature from 
the one shown in Fig. |4] (first state of the first column) to 
one that is highly localized. Looking at the state space 
representation of the eigenstate for V4 = 20 it becomes 
clear that it is a pure normal mode state. More impor- 
tantly, the participation ratios undergo a sharp transition 
around V4 — 12, signaling a significant change in the na- 
ture of the eigenstate. It is intriguing to observe that 
the 156-resonance basis value exhibits a much sharper 
transition as compared to the other basis. 

The various scaled expectation values shown in the in- 
set of Fig. [6] confirm the results based on the participation 
ratios. In particular, an equally sharp transtition is also 
seen in this case at V4 = 12 and hence establishes the 
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IV. CONCLUSIONS 
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FIG. 6. Variation of the participation ratios of the highest 
energy eigenstate in the polyad [7, 14, 1)4] with increasing ex- 
citation of the out-of-plane bending mode 1^4. The zeroth- 
order basis, 156, 526, and 125 single resonance basis values 
are shown as black line, filled circles, grey squares, and tri- 
angles respectively. Inset shows the expectation values (Vise) 
(filled circles), (V526) (grey squares), and (V125) (triangles) for 
the highest energy state as a function of 114. Note the sharp 
transition in both measures around Vi = 12. 



sensitivity of the level-velocities to the changing nature 
of the eigenstate. For V4 > 12 aU the expectation values 
are close to zero and thus predict a pure normal mode 
state. Incidentally, the results here agree well with the 
observation of IVR-protected zeroth-order states around 
V4 — 10 in an earlier study[Hj by Sibert and Gruebele 
(see figure 7 of their paper). However, it is interesting to 
note that the expectation values take on maximum val- 
ues where the participation ratios are minimum. Similar 
observations p"5j were made in an earlier study regarding 
the appearance of new modes in highly excited acetylene. 
Thus, the results of Fig. |6] suggest that there is a bifurca- 
tion in the underlying classical phase space which could 
be responsible for the experimental observation of sharp 
features near the dissociation threshold. Undoubtedly, 
more work is needed to establish the detailed bifurcation 
structure of the classical phase space and corresponding 
fingerprints on the quantum eigenstates. Perhaps, such 
an analysis would resolve the puzzle as to why such sharp 
features can be predicted 8J with an effective Hamiltonian 
without any resonance couplings. 



The central point of this work is to emphasize the fact 
that highly excited eigenstates of multiresonant effective 
Hamiltonians form several sequences, corresponding to 
different classes of dynamics, interspersed amongst each 
other in a complicated fashion. However, the tools used 
in this work are ideal to disentangle the sequences and 
gain insights into the nature of the eigenstates. A re- 
markable aspect that emerges from this and other pre- 
vious studies is worth noting - the expectation values 
shown, for example, in Fig. [5] are quantum objects. Nev- 
ertheless, the observed patterns and sequences in the ex- 
pectation values have a one-to-one correspondence with 
the semiclassical angle space densities of the eigenstates. 
Thus, this establishes the close classical-quantum cor- 
respondence that exists for spectroscopic Hamiltonians. 
At the same time, the observation that state space based 
measures of localization or delocalization of the eigen- 
states are insufficient to identify the sequences suggests 
that phase space is the ideal setting where the true nature 
of the eigenstates are revealed - a conclusion that fortifies 
decades of earlier work [5] on systems with smaller degrees 
of freedom. The semiclassical action-angle representation 
is not all that abstract afterall! 

It is heartening to see that the simple level- velocity ap- 
proach is capable of identifying the eigenstate sequences. 
The present work agrees, to a large extent, with the pre- 
vious analvsis|12| and classification of the eigenstates. In 
particular, the eigenstate classes predicted in the ear- 
lier studies do persist in this study which uses a more 
accurate Hamiltonian. Therefore, the different classes 
of eigenstates and predictions of their influence on the 
IVR dynamics of certain zeroth-order bright states are 
expected to be robust even upon inclusion of other weak 
higher order resonances. The approach adopted in this 
work is easily generalized to systems with any number of 
degrees of freedom with several anharmonic resonances. 
The recently proposed method [16] of lifting the eigen- 
states onto the Arnold web i.e., the network of nonlin- 
ear resonances can lead to further detailed insights into 
the anatomy of certain mixed states (cf. last column of 
Fig.|4|. 

Finally, Fig. |4] presents another crucial question. 
Are there experimental signatures of the eigenstate se- 
quences? Clearly, a state localized in the semiclassical 
angle space representation need not be localized in the 
state space representation. However, since the quantum 
expectation values do pick out the sequences correctly, 
it is of some interest to see if there are any spectral sig- 
natures of such different classes of eigenstate sequences. 
Alternatively, Fig. [6] suggests that new modes, like the 
counter-rotator and local bending modes in acetylene, 
might appear in SCCI2 due to the various bifurcations in 
the underlying classical phase space. Again, there must 
be spectral fingerprints of such new modes. More work 
in this direction needs to be done. 
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